Harmonic mixing in two coupled qubits: quantum synchronization via ac drives 
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Simulating a system of two driven coupled qubits, we show that the time-averaged probability to 
find one driven qubit in its ground or excited state can be controlled by an ac drive in the second 
qubit. Moreover, off-diagonal elements of the density matrix responsible for quantum coherence 
can also be controlled via driving the second qubit, i.e., quantum coherence can be enhanced by 
appropriate choice of the bi-harmonic signal. Such a dynamic synchronization of two differently 
driven qubits has an analogy with harmonic mixing of Brownian particles forced by two signals 
through a substrate. Nevertheless, the quantum synchronization in two qubits occurs due to multi- 
plicative coupling of signals in the qubits rather than via a nonlinear harmonic mixing for a classical 
nano-particle. 



I. INTRODUCTION 

While the desire of building code-cracking quantum 
computers [l|, [2| remains a long standing goal, its pursuit 
has pushed forward the enormous progress achieved in 
quantum mesoscopic physics and quantum nanodevices. 
These efforts have already resulted in the development 
of a new class of mesoscopic devices 0-11] and even new 
types of materials known as quantum metamaterials [7, 

Development of such new devices requires a deep un- 
derstanding of dynamics of single and multi-qubit sys- 
tems driven both by ac signals (e.g., electromagnetic ra- 
diation, or ac voltage and/or currents) and noise. Fol- 
lowing this direction of research, several quantum am- 
plifiers have been recently proposed for one and two- 
qubit [13, systems. The idea d, [3 

was to extend 

the stochastic resonance phenomenon (amplification by 
certain amount of noise) for one or two ac driven qubits 
to enhance quantum coherence. Further analogy |lll4l3l | 
between two-qubit system, a Brownian particle driven 
in periodic substrate and a parametric amplifier has re- 
sulted in a proposal to use two coupled qubits as an am- 
plifier of quantum oscillations. 

This work was motivated by the analogy between 
driven Brownian nanoparticles and a system of coupled 
qubits. It is well-known that an overdamped particle 
driven by two harmonic ac signals through the substrate 
can drift in any desirable direction if frequencies of the 
two drives are commensurate. This effect known as 
harmonic mixing [l3 - |22| has been already observed in 
many systems, including vortices in superconductors (23j . 
nanoparticles driven through a pore |24| . current driven 
Josephson junctions [1^ etc. This suggests an idea that a 
coupled two-qubit system also should exhibit a harmonic 
mixing behaviour, but in contrast to the usual classi- 
cal harmonic mixing for overdamped particles, quantum 
harmonic mixing should be via parametric coupling of 
two drives in the quantum master equation. This ef- 
fect can be used to synchronize quantum oscillations in 
the two qubits and can control the average probability 
for a qubit to stay in either ground or excited state by 
changing the relative phase and/or frequencies of the bi- 



harmonic drive. 

Further, our results can be applied to the case when 
one needs to control qubits, which do not have their own 
control circuitry for the reasons of limiting the decoher- 
ence brought in by such extra elements, or because of ac- 
cessibility (e.g., in case of 2D or 3D qubit arrays (26l [27| 
with control circuitry placed at the boundary viz. sur- 
face of the device). Similar problems arise in the so- 
called indirect quantum tomography (see e.g., Refl28l) or 
in quantum computations with access to a limited num- 
ber of qubits (see e.g., Ref. HI). In all these cases, the 
proposed method of harmonic mixing in qubits allows to 
control the state of the second (not directly accessible 
qubit) by varying the frequency and/or phase of the first 
(accessible) qubit. 



II. MODEL 

In order to describe a two-qubit system we will use a 
Hamiltonian in a spin-representation for each qubit with 
the so called Gx — <Jx coupling: 



2 ^ 



(1) 



where tr^ and cr^ are Pauli matrices corresponding to ei- 
ther the first (j = 1) or the second (j = 2) qubit. The 
tunnelling splitting energies Aj are usually determined 
by the geometry and fabrication details of the specific de- 
vice, while the bias energies Cj can be driven externally. 
For simplicity, we consider two identical qubits; that is, 
we assume Ai = A2 = A. Let us drive the qubits by a 
controlled bi-harmonic drive: 



ei(i) = Ai sin(wit) 
ei(i) = A2 sin(a;2i + 0) 



(2) 



In other words, each qubit is driven by its own signal 
and amplitudes, frequencies and relative phase of these 
two signal can be varied at will. The question arises if 
and under what conditions the second qubit can influence 
the coherence and occupation of the ground and excited 
states of the first one. Therefore, we are interested if 
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the second qubit can be used to control the state of the 
first qubit via dynamic synchronization of their quantum 
oscillations. 

The two-qubit density matrix p can be written as 



P 



E 

,b—0.x.y.z 



"^1 



(3) 



This is a straightforward generalization of the standard 
representation of the single-qubit density matrix expres- 
sion using the Bloch vector; the components liab thus 
constitute what can be called the Bloch tensor. Then, 
the master equation, 



dp 
'dt 



H{t)rp 



(4) 



can be written down directly, using the standard approx- 
imation for the dissipation operator F via the dephasing 
{T ,f,i,T ^2), and relaxation (ri,r2) rates, to characterize 
the intrinsic noise in the system. The master equation 
([3]) can be explicitly written as follows [l^, [HI : 



Also, for simplicity, hereafter we assume that the relax- 
ation rates are the same for both identical qubits, i.e., 
r^i = T^2 = and Fi = F2 = F^, and the tempera- 
ture is low enough, resulting in Zt2 = Zti = 1, where 
Ztj = tanh(Aj/2A:BTj) is the equilibrium value of the 
z-componcnt of the Bloch vector. This set of ordinary 
differential equations is an ideal starting point for numer- 
ical analysis of the dynamics of two driven and dissipative 
qubits as it has been proved before [13, EH- We will simu- 
late these differential equations for two differently driven 
qubits and will study quantum harmonic mixing. 
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FIG. 1: Time-averaged Bloch tensor components {X\) = 
(Hox) (a) and {Zi) = (IIoz) (b) for two coupled qubits driven 
by the two harmonic drives ((2| with parameters: A\ = A2 = 
10, <j) = Q, uj\ = 2\/A2 -I- g'^ and the averaging time inter- 
val 5.6 X 10* < ujtt < 1.4 X 10^ (thus, averaging time T 
was 8.4 X 10*/tJi). Other parameters are the simulation step 
dt = 1.13 X 10~^, the number of simulation steps 5 x lO'"*, 
damping F = 10~^, couphng constant g = 1, the tunnehng 
spUtting energies A = 1. To veryfy our numerical results we 
simulate by using both Euler (open circles) and second-order 
muhidirivative methods (filled circles). The time-averaged 
Bloch tensor element {Xi) responsible for the qubit coher- 
ence shows peaks at bJijiOY = 2/5,4/5,2,4, while the time 
averaged component (Zi) peaks at iiJijiox = 3/5,9/10,3,5. 
Also, pumping of the excited state for incommensurate fre- 
quencies is clearly seen: |(Zi)| increases for L02 > 2bj\. By 
simulating the ten times denser point mesh for frequency ra- 
tios from 0.01 to 1 (see inset in (a)) we have obtained ex- 
tra commensurate resonance frequency ratios but still could 
not resolve width resonances, which is consistent (see, e.g., 
Ref. [21I ) with zero- width harmonic mixing resonances in clas- 
sical nonlinear devices, where the response near resonances 
behave as cos AajT, where T is the observation time, and Ail 
is detuning. In order to check whether we have a similar be- 
haviour for qubit harmonic mixing, we have simulated {Z\) 
and (Xi) (panels (c) and (d)) near resonances bJijiOY — 3 
and 2 with frequency detuning Aoj < 1 /T and observed finite 
resonance width and oscillations near resonances similar to 
classical harmonic mixing. 
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In the limit of zero coupling, 5 = 0, there exists a 
solution of Eqs. ([S]) with no entanglement between the 
qubits. This solution can be written as a direct product 
of two independent density matrices expressed through 
their Bloch vectors: 

The components of the Bloch tensor liab are all zero with 
the exception of 

(Hoa;, Hoy, IIoz) — {Xi, Yi, Zi); {Uxo, Hyo, ^zo) = {X2, Y2, Z2) 

(7) 

which are just the separate qubits' Bloch vector com- 
ponents. The density matrix components Hox = Xi, 
Ilaio = X2, Hoy = Yi, Hyo — Y2, Hoz = Z\, and IIzo — ^2 
can be often directly accessible in experiments. For in- 
stance, for two coupled flux qubits the circulating cur- 
rents in each qubit are proportional to X\ or X2 , respec- 
tively (see e.g., Ref. H), while Z\ and Z2 determine the 
occupation probabilities of the upper (lower) level for the 
first and the second qubit: 



^(lT^,(t)) 



(8) 
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III. SIMULATION RESULTS 



We simulated the set ([5]) by using the standard Eu- 
ler method which has been proved to converge well for 
low-noise drives [lol . [llj and analyzed the time-averaged 
diagonal element of density matrix {X-\) = (^ox) — 
limr-^oo 'HoxdtjT ^ responsible for the mean coherence 
in the first qubit, as well as the time-averaged density 
matrix element {Z\) = (Hoz) = limT^oo Jq ^ozdt/T, re- 
sponsible for the mean occupation of the ground and ex- 
cited states in the first qubit. To verify validity of our 
numerical results we have also used higher-order multi- 
derivative methods to prove the stability of our numerical 
procedure (compare the open circles for Euler methods 
and the filled circles for the second order method in Fig. 
1 and 2). 

As we expected, there is no mean coherence (Xi) w 
for most of frequency ratio uji/lu2 apart from the spe- 
cific commensurate cases (e.g., uji/uj2 — 2/5, 4/5, 2, 4, see 
Fig. la). Such a situation reminds of a usual classical har- 
monic mixing for nanoparticles (see, e.g., plj). however, 
the frequency ratios where peaks occur, are also tune- 
able by changing the absolute value of signal frequency 
in either the first or the second qubit. Indeed, choosing 
the frequency uji to be equal to the inter-level spacing fre- 
quency uji — A'^ + g^—g (Fig. la) or wi — '1\J A? + g"^ 
(Fig. 2a) or even away from the inter-level resonances 
LUi = 2.113(v/A2 + g2 _ g) (Fig. 2b) results in a qual- 
itatively similar peak structure, but showing different 



FIG. 2: Time-averaged Bloch tensor components {Xi) = 
(Iloa;) for two coupled qubits driven by the two harmonic 
drives (Eq.©) with the same parameters as in Fig. 1 and 
driving frequency uji = a/ + — g (a) equal to an energy 
level transition frequency 11] and for cji = 2.113(-\/A^+^— 
(?) (b) which is away from the energy level transition. Simu- 
lations with a ten-time denser point mesh for frequency ratio 
from 0.01 to 1 (see inset in (a)) uncovered some extra com- 
mensurate frequency ratio where peaks in (Xi) occurs. 



sequence of the frequency ratios. Indeed, for uji — 
2.113(\/A2 -t- g"^ — g) (Fig. 2b), several new frequency 
ratios corresponding to the enhancement of qubit coher- 
ence occur at ^2/^1 = 2/5, 4/5, 6/5, 8/5, 2, 12/5, 14/5, 4. 
Moreover, some peaks can even change their signs (com- 
pare peaks at 0^2/^1 = 4/5 in Fig. la, 2a, 2b) indicating 
that both the frequency ratio and the absolute value of 
frequency can be used to tune qubit harmonic mixing. 
Such a behaviour is quite unusual with respect to classi- 
cal harmonic mixing (see e.g., (l6| ) where the frequency 
ratio is defined by nonlinearity of the system. In con- 
trast, the master equation set ([5]) is linear and harmonic 
mixing occurs via a mixture of multiplicative drives as 
in the qubit parametric oscillator. As we have recently 
shown, this results in a quite unusual spectra of Ha,;,, and 
in particular Xi and Zi with many harmonic peaks show- 
ing complex hierarchy. This can explain a non-trivial be- 
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FIG. 3: Dependence of (Xi) and (Zi) on the relative phase 
of two drives of the bi-harmonic signal at different frequency 
ratio bj-ilbj\ — 2 (a), 3(b), 4(c). All other parameters are 
the same as in Fig. 1. For even frequency ratio, where (Xi) 
has a peak (Fig.l), the strong dependence of {X-i){ip) and a 
week dependence of {Z\){ip) occurs, while, for odd ratio of 
0^2/'^!, dependence of {Zi){(p) is clearly seen and {Xi){ip) is 
negligible. The periods of functions (Xi)(ip) and {Zi){(p) are 
controlled by the frequency ratio and are equal to 2-klo\Iloi 
(for a;2 > cJi). 



haviour of harmonic mixing changes when varying ujx or 
ijj2- Note also, that the quantum harmonic mixing occurs 
in both cases: when (i) wi is equal to inter-level spacing 
and (ii) away from this situation. Therefore, there is no 
need to tune the parameters of the external drives to any 
characteristic internal frequency of the two qubit system 
to observe quantum harmonic mixing. 

We have also observed a similar harmonic mixing in 
time-averaged matrix element (Zi) responsible for the 
occupation of the excited and ground states (Fig. lb). 
Interestingly, the peaks in (Zi) occur at different ratios 
of bi-harmonic drive uj^j^i — 3/5,9/10,3,5. However, 
such a behaviour is perfectly consistent with the har- 
monic spectra of (Zi) and (Xi) studied in [il|. Indeed, 
the specrum of X\ contains only odd harmonics, while 
the spectrum of Z\ consist of even harmonics in agree- 
ment with the fact that peaks of (X\) and (Zi) have 



a different parity. Moreover, apart from the peaks at 
the specific commensurate frequencies, the value of |(Zi)| 
gradually increases with the frequency for W2 ^ 2a;i indi- 
cating pumping in the excited state even for incommen- 
surate frequencies. 

Following the analogy with classical harmonic mixing 
[2l| . we expect the dependence of {X]) and (Zi) on the 
relative phase Lp of bi-harmonic drive at commensurate 
frequencies where peaks have been observed. Indeed, we 
obtained such a dependence and (Zi)(<p) shown 

in Fig. 3(a-c) for the same simulation parameters as in 
Fig. 1 and for frequency ratios {jJijhJx = 2 (a), 3 (b), 4 
(c). The well resolved peaks of (Xi) at even frequency 
ratios 0^2/^1 = 2 and 4 exhibit a strong dependence on 
relative phase, while very weak peaks of (Zi) at these 
frequency ratios show almost no dependence on (p. Com- 
paring figures 3(a) and 3(c), we also conclude that peri- 
odicity of the {X\)(lp) changes with increasing frequency 
ratio following the rule: (X\){(p -I- 27ra;i/w4) = {Xi)(Lp). 
Therefore, the number of full oscillations increases with 
frequency ratio W2 /wi . This dependence of the (^-periods 
of (X\i and (Zi) oscillations on the frequency ratio of the 
harmonic drives is analogous to the similar dependence of 
classical harmonic mixing of a Brownian particle driven 
by bi-harmonic drive on nonlinear substrate [2]| . 



IV. CONCLUSIONS 

We have predicted quantum harmonic mixing in a two- 
qubit system driven by a bi-harmonic drive. It manifests 
itself in a set of peaks of time-averaged density matrix 
components responsible for both qubit coherence and 
occupation of ground and excited states of the qubits. 
These peaks can be controlled not only by the ratio of 
frequencies of the two signals but also by tuning frequen- 
cies themselves and by relative phase of the two signals. 
Such a quantum harmonic mixing can be used to manip- 
ulate one driven qubit by applying an additional ac signal 
to the other qubit coupled with the one we have to con- 
trol. Indeed, setting the frequency of the second qubit to 
be three times larger than the one of the first qubit and 
changing the relative phase of signals in these qubits pro- 
duces oscillations of the density matrix element Z\ with 
amplitude of about 0.1 according to Fig. 3b. Therefore, 
by changing the driving signal only in the second qubit 
should allow us to indirectly vary the occupation proba- 
bilities of the upper level in the first qubit between 47 and 
58 percent (see eq. ([H])). A stronger coupling between the 
qubits should allow an even larger amplitude of the indi- 
rect control of the occupation probability. This effect is 
obviously robust to a reasonably strong decoherence and 
energy dissipation in the system. 
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